arXiv:1502.03442v2 [astro-ph.CO] lJul2015 


Mon. Not. R. Astron. Soc. 000, 1—18 (2015) Printed 2 July 2015 


(MN IATeX style file v2.2) 


Disentangling dark sector models using weak lensing statistics 

Carlo R. Benton Metcalf^, Marco Baldi^’^’^, Massimo Meneghetti^’^’^, 

Lauro Moscardini^dd^ Margarita Petkova® 

^ Aix Marseille Universite, CNRS, LAM (Laboratoire d’Astrophysique de Marseille) UMR 7326, 13388 Marseille, Franee 
^ Dipartimento di Fisiea e Astronomia, Alma Mater Studiorum Universitd di Bologna, viale Berti Piehat, 6/2, 40127 Bologna, Raly 
^ INAF - Osservatorio Astronomieo di Bologna, via Ranzani 1, 40127 Bologna, Raly 
^ INFN - Sezione di Bologna, viale Berti Piehat 6/2, 40127 Bologna, Raly 
^ Jet Propulsion Laboratory, 4800 Oak Grove Dr. Pasadena, CA 91109, USA 

® Department of Physies, Ludwig-Maximilians-Universitaet, Seheinerstr. 1, D-81679 Muenehen, Germany 


ABSTRACT 

We perform multi-plane ray-tracing using the GLAMER gravitational lensing code 
within high-resolution light-cones extracted from the CoDECS simulations: a suite of 
cosmological runs featuring a coupling between Dark Energy and Cold Dark Matter. 
We show that the presence of the coupling is evident not only in the redshift evolution 
of the normalisation of the convergence power spectrum, but also in differences in 
non-linear structure formation with respect to ACDM. Using a tomographic approach 
under the assumption of a ACDM cosmology, we demonstrate that weak lensing mea¬ 
surements would result in a ag value that changes with the source redshift if the true 
underlying cosmology is a coupled Dark Energy one. This provides a generic null test 
for these types of models. We also find that different models of coupled Dark Energy 
can show either an enhanced or a suppressed correlation between convergence maps 
with differing source redshifts as compared to ACDM. This would provide a direct way 
to discriminate between different possible realisations of the coupled Dark Energy sce¬ 
nario. Einally, we discuss the impact of the coupling on several lensing observables for 
different source redshifts and angular scales with realistic source redshift distributions 
for current ground-based and future space-based lensing surveys. 

Key words: galaxies: halos - cosmology: theory - dark matter - methods: numerical 
simulations - gravitational lensing: weak 


1 INTRODUCTION 

It is now widely accepted by the scientific community that 
the energy content of our Universe must be largely domi¬ 
nated by some unknown particles and fields beyond the stan¬ 
dard model of particle physics. These are characterised by 
extremely weak interactions with the electromagnetic field, 
and are thereby termed the “dark" components of the Uni¬ 
verse (Suzuki et al. 2012; Planck Collaboration et al. 2013), 
and generically classified as “dark matter" and “dark energy" 
based on their background and clustering properties. While 
dark energy is supposed to source the observed accelerated 
cosmic expansion (Riess et al. 1998; Schmidt et al. 1998; 
Perlmutter et al. 1999; Kowalski et al. 2008) and to have 
at most weak spatial density fluctuations, dark matter con¬ 
stitutes more than 80 percent of the clustering mass in the 
Universe and drives the growth of cosmic structures (Lacey 
& Cole 1993; Tormen 1998; Springel et al. 2005; Giocoli et al. 
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2007) as well as the deflection of light rays as predicted by 
General Relativity (Bartelmann & Schneider 2001; Bartel- 
mann 2010). The latter phenomenon is known as “gravita¬ 
tional lensing": the light traveling from background sources 
(such as distant galaxies) down to the observer is deflected 
by the inhomogeneous intervening matter distribution, caus¬ 
ing a distortion of light bundles and consequently a mod¬ 
ification of the observed galaxy shape (Kaiser & Squires 
1993; Kaiser et al. 1995). Depending on the overall mag¬ 
nitude of the light deflection we can distinguish between 
two main regimes of gravitational lensing: the strong lensing 
(SL) characterising the large distortions generated by single 
highly overdense regions of the Universe (as the central re¬ 
gions of galaxies and galaxy clusters) and the weak lensing 
(WL) that occurs as the integrated effect of light rays trav¬ 
eling through the inhomogeneous cosmic web (Meneghetti 
et al. 2005a; Kneib & Natarajan 2011). In particular, weak 
gravitational lensing represents an important and widely- 
used tool for cosmological investigation and to probe the 
matter density around galaxies and galaxy clusters (Man- 
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delbaum et al. 2006b; Okabe et al. 2010b; Oguri et al. 2012). 
However, as this effect is weak - by definition - it is neces¬ 
sary to average over a large area of sky in order to extract 
a statistically significant lensing signal. 

In the present work we will focus on the prospects for 
using WL as a cosmological probe to distinguish among dif¬ 
ferent models of dark energy. By analysing the distorted 
shape of background galaxies as a function of their redshift, 
WL can in principle constrain the total matter density of 
the Universe, Qm, the linear matter power spectrum nor¬ 
malisation, (78, and the dark energy equation of state icde- 
Furthermore, this could reveal possible signatures of an in¬ 
teraction between the two dark components (Beynon et al. 
2012; Carbone et al. 2013; Pace et al. 2015). In order to 
achieve such discriminating power, however, highly accurate 
and unbiased shear measurements - as the ones expected for 
the next generation of wide surveys - will be required. Since 
gravitational lensing does not depend on the bias between 
the distributions of dark and luminous matter, it represents 
a complementary probe to other constraining observations 
like supernovae, BAO (Baryon Acoustic Oscillations) and 
CMB (Cosmic Microwave Background) measurements (Kil- 
binger et al. 2013; Kilbinger 2014). 

In recent years, the activities of the CFHTLenS collab¬ 
oration have greatly contributed to measuring the cosmic 
shear signal in different patches of the sky with the aim 
to constrain cosmological parameters. From the first explo¬ 
ration by Fu et al. (2008) of weak lensing by large-scale 
structure in the linear regime on a region of 57 square de¬ 
grees, the CFHTLenS work has progressed through a se¬ 
ries of different steps (Benjamin et al. 2013; Kilbinger et al. 
2013) and improvements (Hildebrandt et al. 2012; Heymans 
et al. 2012, 2013) before performing a full weak lensing mea¬ 
surement in three dimensions using a spherical harmonic ap¬ 
proach (Kitching et al. 2014). As demonstrated by the shape 
measurement of galaxies in the COSMOS field observed by 
the Hubble Space Telescope (Amara et al. 2012), in the near 
future, WL space-based measurements (Refregier et al. 2002, 
2004) combined with other independent probes of the large- 
scale matter distribution will be able to discriminate with 
high accuracy among various possible scenarios for the fun¬ 
damental constituents of the dark sector (Harnois-Deraps 
et al. 2014; Kitching et al. 2014b). Finally, WL is also a 
promising tool for identifying clusters (Maturi et al. 2005) 
in blank fields, complementing and possibly driving other 
independent approaches (Bellagamba et al. 2011). 

Ray-tracing through light-cones extracted from numeri¬ 
cal simulations represents the most accurate method to com¬ 
pute WL predictions (Hilbert et al. 2008). In particular, this 
is the case for cosmologies that have not yet been fitted to 
semi-analytical prediction tools for their behaviour in the 
non-linear regime (Lewis et al. 2000; Takahashi et al. 2012). 
The drawback of this approach is that it is computation¬ 
ally demanding as it requires to sample the full cosmolog¬ 
ical parameter space with a series of N-body simulations. 
However, high speed computers and improvements of nu¬ 
merical solver algorithms are progressively simplifying the 
ray-tracing methodology. 

Studying the WL signal signal as a function of source 
redshift is commonly called “lensing tomography” and is 
likely to have great significants for cosmology in the future 
(Schrabback et al. 2010). This is the technique we will in¬ 


vestigate in the present work with the aim of understanding 
whether a tomographic slicing of background sources might 
increase the information extracted from WL observables on 
possible interactions between dark matter and dark energy. 
In particular, the possibility of fully exploiting this method 
to measure specific dark energy signatures does not only rely 
on a high accuracy of the shear measurements, but also on a 
high source density which is necessary to reduce statistical 
errors. 

The Codecs cosmological simulations are hydrody- 
namical simulations but do not include gas cooling, star for¬ 
mation and the feedback mechanism of the baryonic com¬ 
ponent, which may imprint a non negligible signal in the 
cosmic shear power spectrum on angular scales smaller than 
a few arcminutes (Semboloni et al. 2013; Fedeli et al. 2014; 
Harnois-Deraps et al. 2014). However, as we are interested 
in examining relative differences of the coupled dark matter- 
dark energy models with respect to the standard ACDM one 
we will make the common assumption that the relative dif¬ 
ferences between models would be only weakly affected by 
baryonic physics. This assumption is also reinforced by the 
fact that the baryonic physics has a non-negligible impact 
for scale k > 3500/i/Mpc. 

Pace et al. (2015) have already performed ray-tracing 
simulations in different coupled dark matter-dark energy 
models extracted from the CoDECS suite, constructing maps 
for different lensing quantities starting from the lensing po¬ 
tential. Their analysis is performed on maps with a resolu¬ 
tion of around 20 arcsec and considering sources located at 
redshift Zg = 1. They found that the most significant differ¬ 
ences from the standard ACDM model are due to differences 
in the growth of the perturbations and to the effective fric¬ 
tion term in non-linear dynamics. The most extreme realisa¬ 
tion of coupled Dark Energy expectedly showed the largest 
difference from ACDM of about 40% in the power spectrum, 
as found also by Carbone et al. (2013). In this paper we 
extend and complement their analyses focusing the atten¬ 
tion on the small scales regime thanks to high-resolution 
ray-tracing simulations performed with the GLAMER code. 
In addition we construct ray-tracing simulations considering 
different source redshift distributions, with the aim of to- 
mographically evaluate the difference between coupled dark 
energy models and the standard ACDM one. 

The structure of this paper is organized as follows. In 
Section 2 we summarise the numerical simulations that are 
used to construct the light-cones and in Section 3 we de¬ 
scribe the methodology we adopted to perform multi-plane 
ray-tracing. Results and lensing statistical analyses are pre¬ 
sented in Section 4. Einally, in Section 5 we discuss our con¬ 
clusions and summarise our main results. 


2 COSMOLOGICAL SIMULATION 

Eor our investigation we rely on the publicly available 
CoDECS simulations (Baldi 2012c) that represent the largest 
suite of cosmological and hydrodynamical simulations of 
coupled Dark Energy (cDE) models to date. The simulations 
have been performed by means of a modified version devel¬ 
oped by Baldi et al. (2010), of the widely used TreePM/SPH 
N-body code GADGET (Springel 2005), and self-consistently 
include all the effects associated with the non-minimal inter- 
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Table 1. The list of cosmological models considered in the present work and their specific parameters. All the models have the same 
amplitude of scalar perturbations at zcmb ~ 1100, but have different values of ag dit z = 0. In short, a is a parameter in the inflation 
potential as shown, /3(0) is the coefficient of the coupling term with dark matter density and W(i){z = 0) is the effective equation of state 
parameter (p/p). See Baldi (2012c) for details. 


Model 

Potential 
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W(f,{z = 0) 

asiz = 0) 
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-1.0 
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v(<j)) = 
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Figure 2. An illustration of the construction of the light-cone up to redshift 2 ; = 1.4 with an aperture of 5 X 5 square degrees for the 
ACDM model. The light-cone passes through three simulation boxes, since the comoving distance corresponding to 2 ; = 1.4 is 3 Gpc/h. 
As indicated in the upper part of the flgure we have 12 snaphosts available up to this considered source redshift from which we construct 
14 lens planes. In the square that contains a single simulation box we stack together slices from difference snapshots of the box. The 
snapshot id numbers are shown on the top of each rectangle representing the slices. The different colors used to shade the squares indicate 
that the simulation boxes have been randomized as described in Roncarelli et al. (2007). The pink shaded triangle shows the region of 
the simulation snapshots used to construct the planes within the light-cones. 


action between a DE scalar field cj) and CDM particles. The 
Codecs suite includes several different possible combina¬ 
tions of the scalar DE potential - the exponential (Lucchin 
& Matarrese 1985; Wetterich 1988) or the SUGRA (Brax & 
Martin 1999) potentials for example - and of the coupling 
function which can be either constant or exponential in the 
scalar field (see e.g. Baldi et al. 2011). In the present work, 
we will consider four models (ACDM, EXP003, EXP008e3, 
and SUGRA003) with different combinations of these free 
functions that are summarised in Table 1. Eor more details 
on the models we refer to Baldi (2012c). 

Such variety of scenarios is reflected in the diversity 
of effects that they determine on both the background ex¬ 
pansion history and the linear and non-linear evolution of 
perturbations. More specifically, the dynamical evolution of 
the DE scalar field through a matter dominated scaling solu¬ 
tion alters the cosmic expansion history with respect to the 


standard ACDM cosmology in a way that substantially de¬ 
pends on the potential and coupling functions adopted. This 
is shown in Eig. 1, where we show the comoving distance as 
a function of redshift for the four models. The bottom panel 
displays the relative differences of the comoving distances at 
a given redshift between the cDE models and the reference 
ACDM. Erom the figure we can notice that while locally 
the comoving distances are consistent with each other (due 
to the same normalisation of Hubble parameter, i^o, for all 
the cosmologies), at higher redshifts the cDE models have 
smaller volumes than ACDM, with a maximum effect for the 
SUGRA003 cosmology which deviates low from ACDM by 
about 5 percent already at z = 2. 

At the level of linear density perturbations, the mod¬ 
els all predict an enhanced growth rate with respect to 
ACDM at high redshifts. However, while the EXP003 and 
the EXP008e3 models show an enhanced growth also at 
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Figure 1. Comoving distances as a function of redshifts for the 
four considered models of the CoDECS suite. The bottom panel 
shows the relative difference of the models with respect to ACDM. 


low redshifts, thereby resulting in a larger value of erg at 
z = 0, the SUGRA003 cosmology is characterized by a 
slower growth as compared to ACDM for z < 7 resulting 
in a comparable value of erg (see the last column of Table 1). 
The non-linear effects of these models have been studied in 
several publications based on the outcomes of the CoDECS 
simulations and range from the impact of cDE on the abun¬ 
dance and structural properties of halos (Baldi 2012a; Baldi 
et al. 2011; Cui et al. 2012; Giocoli et al. 2013), on the sta¬ 
tistical properties of the large-scale structures distribution 
(Marulli et al. 2012; Moresco et al. 2014), on the proper¬ 
ties of the Inter-Galactic Medium at high redshifts (Baldi 
& Viel 2010), and on weak leasing statistics (Beynon et al. 
2012; Carbone et al. 2013; Pace et al. 2015). In the present 
work, we aim at extending the latter analysis by investigat¬ 
ing whether a tomographic slicing of the background sources 
within the light-cone of a 25 square degrees field of view 
might provide additional information to observationally dis¬ 
tinguish the cDE models from the standard ACDM cosmol¬ 
ogy and possibly the different cDE models from each other. 

For our analysis we will make use of the L-CoDECS se¬ 
ries consisting of a periodic cosmological box of 1 Gpc//i 
aside filled with 2 x 1024^ particles evolved through colli¬ 
sionless dynamics from z = 99 to z = 0. All the models 
share the same initial conditions at the redshift of the CMB 
^CMB ~ 1100 and have the same cosmological parameters 
at z = 0 consistent with the WMAP7 cosmological results 
(Komatsu et al. 2011), namely Qcom = 0.226, Qde = 0.729, 
h = 0.703, A, = 2.42 x 10“^ fib = 0.0451 and n, = 0.966. 
The mass resolution is Mcdm{z = 0) = 5.84 x 10^° Mq/H 
for CDM particles and Mb = 1.17 x 10^° Mq/H for the 
(collisionless) baryonic particles (see Baldi 2012c, for a de¬ 
tailed discussion), while the gravitational softening was set 
to Cg = 20 kpc//i. 


3 LENSING PIPELINE 

In the following sections we will present the procedure we 
followed in constructing the lens planes from the cosmologi¬ 
cal simulations, assembling them into light-cones and tracing 
the paths of light through them. 


3.1 Constructing the light-cone: MapSim 

Our code for extracting the particles from the simulation’s 
snapshot files and assembling them into a light-cones is 
called MapSim. The steps MapSim goes through in con¬ 
structing a light-cone can be summarized as follows: 

• Read in an input parameter file that contains informa¬ 
tion about the desired field of view, highest source redshift 
(in this case taken to be = 4) and locations of snapshot 
files. The number of lens planes required is decided ahead 
of time in order to avoid gaps in the constructed light-cone. 
The choice Zs = 4 has been made to better understand where 
the dynamical evolution of the DE scalar field and the en¬ 
hanced growth rate - in the different models - start to leave a 
mark in the weak leasing observables. Notice that at low red¬ 
shifts the models show a different behaviour: while EXP003 
and EXP008e3 continue their enhanced growth, SUGRA003 
does not. 

• Read in each snapshot file going from the present time 
to higher redshift snapshots while extracting only the parti¬ 
cle positions within the desired field of view. Only a single 
snapshot is in memory at any time. 

• Selection and randomization of each snapshot is done 
as in Roncarelli et al. (2007). If the light-cone reaches the 
border of a simulation box before it has reached a redshift 
range where the next snapshot will be used, the box is re¬ 
randomized and the light-cone extended through it again. 

• The lensing planes are built by mapping the particle 
positions to the nearest pre-determined plane, maintaining 
angular positions, and then pixelizing the surface density 
using the triangular shaped cloud (TSG) method (Hockney 
& Eastwood 1988). The grid pixels are chosen to have the 
same angular size on ah planes. The lens planes have been 
constructed each time a piece of simulation is taken from 
the stored particle snapshots; their number and frequency 
depend on the number of snaphosts stored while running the 
simulation. 

In Fig. 2 we show an illustration of the construction of the 
light-cone up to redshift = 1.4, for the AGDM model, pil¬ 
ing one on top of the others the different portions of the sim¬ 
ulation snapshots. The differently colored squares are differ¬ 
ent realizations using the same randomization method. The 
rectangles within the squares represents the portion of the 
simulation snapshot from which the particle density distri¬ 
butions have been taken. The number on the top of them 
are the id numbers of the corresponding simulation snap¬ 
shots. The vertical line in the middle of each rectangle indi¬ 
cates the planes onto which the particles are projected for 
ray-shooting. We remind the reader that to save disk space 
not all 92 simulation snapshots have been stored in running 
the simulation. Up to redshift Zg = 4 we saved 18 snapshots 
which are sufficient to consistently model the matter density 
distribution in our light cones. 
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Table 2. Comoving distances in Mpc/h for three different source 
redshifts in the considered cosmological models 


model 

Zs = 0.5 

= 1.4 

Zs — 4 

ACDM 

1327.27 

3000 

5179.52 

EXP003 

1324.56 

2999.75 

5135.05 

EXP008e3 

1320.86 

2981.97 

5091.15 

SUGRA003 

1297.07 

2892.85 

4962.47 


The maps were constructed with a 5 x 5 sq. deg. field 
of view and an angular resolution of 8.8 arcsec. For each 
model we have constructed twenty-five independent realiza¬ 
tions, being careful that in each realization the same field of 
view is selected in the different cosmologies. However, con¬ 
sidering the different comoving distance-redshift evolution 
through the simulation snapshots we have built 22 planes up 
to redshift z = 4 for the ACDM, EXP003 and EXP008e3, 
while we have 21 planes for the SUGRA003 model. In Ta¬ 
ble 2 we point out the comoving distances corresponding to 
redshift Zs = 0.5, 1.4 and 4 for the four models. 



3.2 Ray-Tracing through the planes 


Once the lens planes are created as described in the previ¬ 
ous section the lensing calculation itself is done using the 
GLAMER lensing code (Metcalf & Petkova 2013; Petkova 
et al. 2013). The multiplane ray-tracing method is described 
in detail in Petkova et al. (2013) so we will only outline the 
procedure here. 

A few definitions are required. If the angular position on 
the sky is 6 and the position on the source plane expressed 
as an angle (the unlensed position) is then a distortion 
matrix A can be defined as 

A= ^ - 

de~ 

The traditional decomposition of this matrix is shown, where 
K is called the convergence and 7 represents the shear. The 
torsion, 73 , represents a rotation which can occur when there 
are multiple deflection planes. It is of order the shear squared 
(see Petkova et al. 2013) and according to our numerical 
calculations, and those of Becker (2013), it is quite small, 
but it will be retained here for completeness. 

When there is a single lens plane, the convergence can 
be expressed as a dimensionless surface density. 
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Xcrit — 


Di 


( 2 ) 


( 3 ) 


47rG DsDis 

is called the critical density, c is the speed of light, G is New¬ 
ton’s constant and Di Ds and Dis are the angular diameter 
distances between observer-lens, observer-source and source- 
lens, respectively. In general, with multiple lens planes, this 
is not the case however. 

The deflection caused by a lens plane, a, is related to the 
surface density on the plane, X(x), through the differential 
equations 

V^^(x) = ^^S(x) , a(x) = V4>W). (4) 


i -0.5 r 
o 

2 ? 102 
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10 ^ 
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Figure 4. Convergence power spectrum of the four considered 
Codecs cosmologies for three different source redshifts: Zs = 
0.5, 1.4 and 4. The shaded grey region encloses the standard devi¬ 
ation of the mean associated to the ACDM model on the different 
5x5 degree light-cone realisations. In the top panel the solid and 
dot-dashed curves show the linear and non-linear CAMB predic¬ 
tions for the ACDM cosmology respectively. For the non-linear 
power spectrum prediction we adopt the extended version of the 
Halofit Model (Smith et al. 2003) from Takahashi et al. (2012). 
In the bottom panels we present the relative residuals of the con¬ 
vergence power spectra with respect to the ACDM cosmology of 
the three cosmological models - for each of the considered source 
redshift cases - featuring a direct interaction between the Dark 
Energy and Dark Matter. 


where the derivatives are with respect to the position on 
the lens plane. These equations are solved on each source 
plane by performing a Discrete Fourier Transform (DFT) 
on the density map, multiplying by the appropriate factors 
and then transforming back to get a deflection map with the 
same resolution as the density map. With the same DFT 
method the shear caused by each plane is simultaneously 
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Figure 3. Convergence maps of a light-cone with aperture 100 arcmin on the x-side, for the four considered models - in each panel top 
left ACDM, top right EXP003, bottom left EXPOOSeS and bottom right SUGRA003. The three panels refer to different cases: sources 
located at redshifts Zs = 0.5, 1.4 and 4, from top to bottom, respectively. In each map the sticks show the direction of the corresponding 
field. 
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calculated. Since the rays are propagated between planes 
using the standard distances in a Robertson-Walker metric 
which assumes a uniform distribution of matter the addition 
of matter on each of the planes will, in a sense, over-count 
the mass in the universe. Without correcting for this the 
average convergence from the planes will be positive and will 
cause the average distance for a fixed redshift to be smaller 
than it should be. To compensate for the implicit density 
between the planes the ensemble average density on each 
plane is subtracted. Each plane then has zero convergence 
on average and the average redshift-distance relation is as it 
would be in a perfectly homogeneous universe. 

After the deflection and shear maps on each plane 
are calculated the light-rays are traced from the observers 
through the lens planes out to the desired source redshift. 
The shear and convergence are also propagated through the 
planes as detailed in Petkova et al. (2013). GLAMER per¬ 
forms a complete ray-tracing calculation that takes into ac¬ 
count non-linear coupling terms between the planes as well 
as correlations between the deflection and the shear. No 
weak leasing assumption is made at this stage. The rays 
are shot in a grid pattern with the same resolution as the 
mass maps: 5 degrees resolved with 2048 pixels on a side. 

In Eig. 3 we show the convergence maps of the same 
light-cone realisation extracted from the different models - 
top left ACDM, top right EXP003, bottom left EXP008e3 
and bottom right SUGRA003. We show the maps for sources 
located at three fixed redshifts: Zs = 0.5, 1.4 and 4 from top 
to bottom, respectively. The sticks in each panel indicate the 
directions of the corresponding shear field. As discussed in 
Gui et al. ( 2012 ) and Giocoli et al. (2013), we immediately 
notice that the density distribution differs in cDE models 
from AGDM due to a difference in the growth as a function 
of redshifts. The top panels show the presence in the field 
of view of a cluster, at redshift z < 0.5, which, while it ap¬ 
pears “assembled” and with a single peak in the AGDM and 
the SUGRA003 models, it is less evolved, showing multiple 
components, in EXP003 and EXP008e3. The intermediate 
and the high redshift maps also exhibit differences mainly 
because of differences in the evolution of the power spectrum 
normalisation and non-linear structure formation. 


3.3 Analytic methods 

An approximation is commonly used to calculate the con¬ 
vergence power spectrum that avoids the complications dis¬ 
cussed in the previous two sections. The convergence can 
be calculated by adding up the single plane convergences, 
equation (4), along an unperturbed light-ray, the Born ap¬ 
proximation. This results in the expression 


k{Zs,0) 


3HiQm D{z)D{z,zL 

2c2 io D{zs)aizy 


5(D{z)e,z) , 
( 5 ) 


where w is the radial comoving distance, a = (l+z)“^ is the 
scale factor and (5(x) = (p(x) — p)/p is the density contrast. 
The angular power spectrum of n in the small angle limit is 
then found to be 


PW) = 


4c4 



/ Diz,Zs) 
\D(zs)a{z) 


2 

Ps 


. D{z) 


(6) 


(Kaiser 1992). Either an analytic model for the density 
power spectrum, Ps{k,z), or a power spectrum taken di¬ 
rectly from a simulation can be inserted into eq.( 6 ). 

In the following sections we will contrast our direct cal¬ 
culations with some analytic models using this approxima¬ 
tion. Eor some applications this approximation is adequate, 
but in new situations and for particular statistics it needs 
to be checked against simulations. 


4 RESULTS 

4.1 Convergence Power Spectrum 

In the real (angular) space, a direct measurement of weak 
leasing is the two-point shear correlation functions and 
that can be obtained by averaging over galaxy pairs with 
angular separations \0i — 0j\ within a bin 6: 

Si(f') 7 (*) 

where the measured galaxy ellipticity measurements et and 
Cx are the tangential and cross components with respect 
to the line connecting the pair, respectively. The weights w 
are obtained from the galaxy shape measurement pipeline 
(Bernstein & Jarvis 2002; Bridle et al. 2010; Kacprzak et al. 
2012; Miller et al. 2013). These two-point shear correla¬ 
tion functions can be calculated from the convergence power 
spectrum by the relation: 

^+/_ = ^ d/ lPW)Jo/4m , (8) 

where Jo and J 4 are the Bessel functions and we have set the 
B-mode power spectrum to zero because leasing generates 
only E-modes in the weak leasing limit. 

In Eig. 4 we show the convergence power spectrum up 
to / ~ 10 ^ measured in the four cosmologies, considering 
sources at three different redshifts: Zg = 0.5, 1.4 and 4. The 
curve referring to each model represents the average over 
all the light-cones. Eor the AGDM case the shaded region 
encloses the standard deviation of the mean associated to 
the different realizations. The Imin considered corresponds 
to the minimum resolvable in the assumed field of view. 
In the same figure we also show the predictions obtained 
by inserting into equation ( 6 ) analytic models for the lin¬ 
ear and the non-linear power spectrum, namely the Halofit 
Model (Smith et al. 2003) from Takahashi et al. (2012) im¬ 
plemented in GAMB (Lewis et al. 2000). As shown in Pace 
et al. (2015) - where they consider the matter density dis¬ 
tribution only up to redshift Zg = 1 - the EXP003 model 
has much more power than the AGDM with the largest dif¬ 
ference being on small scales and increasing with the source 
redshift. Interestingly, for the SUGRA003 model, the power 
spectrum is below AGDM for low multipoles and above it 
for high multipoles. This is a result of the fifth force term 
in the dark matter-dark energy coupling which drives rapid 
structure formation at high redshift, but slows it down at 
lower redshift. The SUGRA003 and the AGDM power spec¬ 
tra intersect at about I ~ 10 ^, slightly decreasing with the 
source redshifts. Erom Eig. 4 we also see that the EXP008e3 
model has about 15-25 percent more power than AGDM 
at all scales and source redshifts, a signature of its growth 
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Figure 5. Left panel: ratio between the convergence power spectrum computed at different source redshifts - as indicated in the labels 
- and the one at 2:5 = 0.5 for each cosmological model. Right panel: rescaled cross power spectra between redshift 2 : = 0.5 and Z 2 - as 
indicated in the labels - for the different models, computed from the lensing maps. All the cross spectra have been rescaled with respect 
to the ACDM prediction. This is why the correlation parameter on the right can be greater than one. Different line styles and colors are 
as in Fig. 4. 


rate being maintained over the whole redshift range. As ob¬ 
tained by the tomographic analysis performed by Pace et al. 
(2015), where only the ratio of the power spectra for sources 
at redshift Zs — ‘I and Zs = 1 has been considered, we con¬ 
firm that the small differences they find at very high I in 
the coupled models, and particularly in the bouncing model 
SUGRA003, are actually present and are much clearer in 
our high-resolution ray-traced maps. The same behaviour of 
the power spectra we find has also been noticed and dis¬ 
cussed by Carbone et al. (2013) using ray-tracing technique 
to compute the CMB lensing maps; however also Carbone 
et al. (2013) are not able resolve the small scale features of 
the coupled DE-DM models because the angular resolution 
on their maps is more than a factor of ten lower than ours. 

Clearly weak lensing tomography, i.e. the study of the 
weak lensing signal as a function of the source redshift, 
can be an important tool for studying these differences 
in the evolution of structure with redshift that occur in 
different coupled Dark Energy models. Eig. 5 shows this 
more clearly. In the left panel we display the ratio between 
the convergence power spectrum computed at five different 
source redshifts and the one computed for sources located 
at 2^5 = 0.5; line styles and colors are the same as in Eig. 4. 
Eor EXP003 and EXP008e3 the ratios tend to lie above the 
ACDM one, but for SUCRA003 it stays below. This is a 
result of structures in the SUCRA003 model evolving less 
rapidly at late times while in the EXP003 and EXP008e3 
models they evolve more rapidly. The difference between the 
models is most evident at small scales as a consequence of 
non-linear structure formation. The right panel of Eig. 5 
show the rescaled cross-spectra between the convergence at 
different redshifts within the same light-cone. High correla¬ 
tion indicates that the lensing is being caused by the same 
objects. Interestingly, SUGRA003 has more correlations at 


small scales and less at large scales, while the other mod¬ 
els have the opposite trend, due to an enhancement of their 
growth rate at high redshifts and to a depletion at low red¬ 
shifts with respect to the ACDM one. We remind the reader 
that the growth rate in EXP003 and EXP008e3 is always en¬ 
hanced with respect to the standard model causing a higher 
as at the present time. 

The pixel resolution of our maps (8.8 arcsec in all cases) 
will have some impact on the accuracy of the lensing statis¬ 
tics calculated. In particular the pixelation tends to smooth 
out peaks and reduce the power on small scales (Takahashi 
et al. 2011; Pace et al. 2015). To investigate how our pix¬ 
elation might be affecting the results we constructed the 
pixel Probability Distribution Punction (PDE) of the con¬ 
vergence map as shown in Eig. 6 for one realisation of the 
ACDM model with sources at Zs = 1.4. The different line 
styles and color histograms show the PDE when the orig¬ 
inal map has been pixel-degraded by a factor of 2, 4 and 
8 - see the figure caption for more details. It can be seen 
that increasing the pixel size reduces the number of very 
high and very low k pixels as expected. The impact this has 
on the convergence power spectrum is shown in the right 
hand panel of Eig. 6, where the power spectra for different 
resolutions are plotted. In the bottom right hand panel are 
the relative residuals with respect to the original, highest 
resolution map. These figures show that reducing the map 
resolution by a factor of two causes small differences of only 
a few percent in the power spectrum relative to the original 
map for / ^ 5 x 10^. The discrepancies become larger for 
the maps degraded by a factor of four and eight. We con¬ 
clude that our calculations are not affected by pixelation for 
/ ^ 5 X 10^ above few percent level. 

It has been shown in this section that the evolution 
of the density power spectrum in different coupled dark en- 
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Figure 6. Left panel: histograms of the convergence values in one realization of the ACDM model, considering sources at 2:5 = 1.4. The 
solid black histogram shows the case for the map resolution used throughout this paper ( 8.8 arcsec pixels). The short dashed green, 
long dashed magenta and dot-dashed blue histograms show the PDF of the maps degraded by a factor of 2, 4 and 8 , respectively. Right 
top panel: convergence power spectra of these maps. Right bottom panel: the ratios of the power spectra with different resolutions with 
respect to the highest resolution one. 
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Figure 7. Left panel: ag as a function of the source redshift, obtained by fitting the corresponding convergence power spectrum for each 
cosmological model, fixing the total mass density parameter Dm to its true value in the simulation. The colored points show the value of 
as at the present time interpolated from linear theory: blue triangle (EXP003) , orange diamond (EXP008e3), red cross (SUGRA003) 
and black dot (ACDM). Right top panel: convergence power spectra for sources at redshift 2 ;^ = 4 as measured from the light-cones in the 
four cosmological models. As in previous figures solid, dotted, dashed and dot-dashed lines correspond to ACDM, EXP003, EXP008e3 
and SUGRA003, respectively. The corresponding solid colored curves represent the best ACDM power spectra obtained fixing the total 
matter content and varying ag- Bottom right panel: the ratio between the convergence power spectrum with the best ag and the one 
obtained using ag(z = 0), i.e. the interpolated value using linear theory (Baldi 2012c) is shown for the three cDE models. 
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ergy models will cause significant changes to the convergence 
power spectrum as a function of source redshift. These differ¬ 
ences from ACDM can be positive or negative at a particular 
source redshift and they can be scale dependent. In the fol¬ 
lowing sections we will investigate some other statistics that 
might be observationally more practical in terms of direct 
measurements and noise estimations. 

4 . 1.1 Impact on the measured normalisation of the power 
spectrum for sources at different redshifts 

As noted previously, the coupled DM-DE models affect both 
the power spectrum normalisation and the small scale non¬ 
linear behaviour of structure formation when compared with 
the ACDM model. One way to see this is through the im¬ 
pact on measurements of the power spectrum normalisation. 
The power spectrum normalisation will be quantified by the 
standard ag parameter which is the variance of the mass 
overdensity density within a sphere of radius S h~^ Mpc 
at redshift zero with only the linear evolution of the power 
spectrum taken into account. When ag is measured using 
data at redshifts larger than zero a correction factor that is 
cosmology dependent must be applied to translate the nor¬ 
malisation to z = 0. If the underlying cosmology assumed 
in doing this procedure is the correct one then ag will not 
depend on the source redshift. In this section we investigate 
what would happen if a ACDM cosmology is assumed while 
the true cosmology corresponds to one of our coupled dark 
energy models. This exercise will help us gain some insight 
on when during cosmic history these models leave most of 
their imprints on the lensing power spectrum. 

The left hand panel of Fig. 7 shows the recovered ag 
values as a function of the different source redshifts for the 
four cosmological models while assuming ACDM. We mea¬ 
sure ag integrating the non-linear power spectrum (see e.g. 
Lewis et al. 2000; Takahashi et al. 2012) up to the differ¬ 
ent considered source redshifts assuming the ACDM model; 
in this way the ag value is directly related to the measure¬ 
ment of the growth factor associated to a ACDM Universe. 
Line styles and colors are the same as in Fig. 4, the shaded 
light grey and grey regions indicate 5 and 10 percent accu¬ 
racies in the measurement of ag for the ACDM case. Also 
shown are the true ag values - interpolated using linear the¬ 
ory (Baldi 2012c) - for each model. In the case of ACDM 
the correct value is recovered, but in the other models incor¬ 
rect ag values are recovered and they change as a function of 
source redshift. Such behaviour would therefore signal a fail¬ 
ure of the underlying ACDM assumption. In EXP003 and 
EXP008e3, ag is underestimated, while the opposite is true 
in SUGRA003 for some of the source redshift range. These 
trends are due to the different evolution of the growth factor 
and of the Hubble function in the coupled DM-DE models - 
see Fig. 1 and Fig. 2 in Baldi (2012c). It is interesting to note 
that for the three coupled models we never obtain the value 
of ag interpolated from linear theory at the present time. 
This means that when fitting the convergence power spec¬ 
trum of the coupled DM-DE models both projection effects 
and non-linearities contribute to the measured ag. 

In the right hand panel of Fig. 7 we show the mea¬ 
sured convergence power spectrum in the different light- 
cone simulations for sources sot Zs = 4 - solid, dotted, 
dashed and dot-dashed refer to ACDM, EXP003, EXP008e3 


and SUGRA003, respectively. For the three coupled DM-DE 
models we display best fitting theoretical convergence power 
spectra obtained by assuming we live in a ACDM universe. 
In the bottom panel we present the relative deviation of the 
theoretical convergence power spectrum computed with the 
best ag and the one linearly interpolated from theory at 
z = 0. While for SUGRA003 model the relative deviation is 
of the order of 1 — 2%, for both EXP003 and EXP008e3 it 
manifests larger values. 

These results indicate that, if sufficient accuracy is at¬ 
tained, measuring ag from shear maps as a function of red¬ 
shift would lead to inconsistencies were one of these coupled 
dark energy models the correct one. To observe these ef¬ 
fects will require large amounts of very accurate data of the 
kind that will be provided by future wide surveys, like the 
Euclid space mission (Laureijs et al. 2011). An additional 
complication not taken into account here is that weak lens¬ 
ing studies are generally sensitive to a combination of the 
power spectrum normalisation, ag, and the total density of 
matter, Qrn, in the combination ag (Qrn/0.25)^ where a is 
dependent on the weak lensing statistic used: a = 0.46, 0.53 
and 0.64, for the shear two-point correlation function, the 
shear in a top-hat or aperture mass, respectively (Fu et al. 
2008; Kilbinger et al. 2013, see also a more extended discus¬ 
sion in the next section). However, is likely to be well 
constrained by other observations such as the CMB. 

4.2 Other weak lensing statistics 

Other statistics of the convergence field besides the power 
spectrum may help to probe the non-Gaussian nature of 
the probability distribution function of gravitational lensing 
observables. In particular, statistics like the variance and 
the skewness in a top-hat or compensated filter represent 
interesting tools to constrain cosmological parameters and 
the dark energy evolution as a function of redshift (Fu et al. 
2008; Kilbinger et al. 2013; Kitching et al. 2014). In this 
context, multi-plane lensing simulations are important tools 
to compute the predicted variance and skewness of the shear 
or convergence field, both for standard and non-standard 
models. In particular, the skewness in aperture is somehow 
independent of the normalization of the power spectrum and 
represents a strong indicator of the evolution of Urn as a 
function of redshift (Schneider et al. 1998). 

The interpretation and modelling of the lensing signals 
due to large scale structures requires a precise understand¬ 
ing of the high-order statistics of the lensing field (Sato 
et al. 2009). Several cosmological analyses have been done 
by studying the mass aperture variance and the skewness 
as a function of the smoothing angle 0, adopting a top-hat 
or a compensated aperture filter. These measured quanti¬ 
ties have the advantage that they can be directly compared 
with the theoretical predictions extracted from the conver¬ 
gence power spectrum. However, a good model of the signal 
on small scales, where the non-linear effects start to domi¬ 
nate, is of fundamental importance to interpreting observa¬ 
tional data. Since the top-hat shear rms between different 
CoDECS models, but only at Zs = 1, has been studied by 
Pace et al. (2015), we will present and discuss in this sec¬ 
tion the aperture-mass dispersion and the associated skew¬ 
ness for three different source redshifts. 

The variance of the filtered shear field as a function of 
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Figure 8. Top panels: variance of the convergence field smoothed with a compensated aperture filter as a function of the smoothing 
scale for sources at three different redshifts: 2:5 = 0.5, 1.4 and 4, from left to right. Bottom panels: relative residuals as a function of 
the smoothing scale with respect to the measurements in the ACDM model. The grey shaded region represents the rms of the variance 
computed on the different realizations of the ACDM light-cone. In the right panel, for sources at Zs = 4, we also show predictions from 
linear and non-linear theory. For the other source redshifts the agreement is analogous to this one. 


the filter size contains the same information as the power 
spectrum, but the filter can be made to have compact sup¬ 
port and additional practical advantages (Sato & Nishimichi 
2013). They can also have a different sensitivity to the 
cosmological model and be particularly dependent on non¬ 
linear structure formation. In addition, the skewness and 
higher-order statistics can be easily defined and interpreted 
for the filtered shear field. We will investigate a particular 
choice of filter called the compensated aperture filter: 




6 

7r0‘^ 




with support 'd = [0,^]. In Fourier space it is 


Wap = 


e 


(9) 


( 10 ) 


the ACDM ones. We notice that the variance of the mea¬ 
surement in the ACDM model decreases as a function of 
the source redshift, as expected from the evolution of the 
clustering of the dark matter in the Universe. While at low 
redshifts the Universe is more clustered and so we may trace 
rays through clusters and voids - enhancing the variance 
of the measurements in different realisations, at high red¬ 
shifts the Universe is more homogenous and so the variance 
through different realisations is expected to be smaller, as 
we confirmed in Fig. 8. 

To demonstrate how much of this signal depends on 
non-linear structure formation we show {M\^) calculated 
with only the linear power spectrum in the right hand panel 
of Fig. 8 for Zs = 4 in the ACDM model. It can be seen that 
non-linear structure formation has a dominant role below 
^ ^ 20 arcmin. 


as Fourier counterpart (Harnois-Deraps et al. 2012). It has 
the attractive property of being well-localized in Fourier 
space near wavenumber I ~ 5/0. We apply this filter in 
Fourier space with zero-padding to reduce boundary effects. 
Since the average value of the convergence on each plane 
is zero, the variance across the map, typically indicated as 
{M^p), can be computed from the convergence power spec¬ 
trum performing the following integral: 

{k^)Ap = (Mlp) = f dllPW)W%p{W) ■ (11) 

(271-) J 

Similarly, it is possible to define the large scale structure 
noise (Hoekstra 2003) adopting a different compensated fil¬ 
ter (see Appendix A for more discussion about this). 

In Fig. 8, we show the variance of the convergence field 
as a function of the smoothing scale 0 with this filter. The 
different line styles and colors refer to the four considered 
cosmological models and the shaded grey region encloses 
the rms of the measurement performed in the different re¬ 
alizations of the ACDM light-cones. The measurements are 
presented for three source redshifts: Zg = 0.5, 1.4 and 4. 
The bottom panels show the relative residuals of the mea¬ 
surements in the different coupled models with respect to 


The agreement between non-linear theoretical predic¬ 
tions and our simulations for ACDM are quite good aside 
from limitations in the simulations at small scales due to 
numerical mass and force resolution limitations. As with 
the power spectrum, the deviations from ACDM in Fig. 8 
are positive for EXP003 and EXP008e3 (by about 50 and 
20 percent - as for the convergence power spectrum). Eor 
SUGRA003, {M\p(0)) is larger than for ACDM at small 0 
and becomes smaller at large scales mirroring the behaviour 
in Eig. 4. This trend is a consequence of a slightly lower as 
governing the large 0 behaviour while for small scales the 
boost in variance comes from differences in the structure of 
halos: in SUGRA003 haloes are more concentrated and have 
more substructures mainly at low redshifts (as was shown 
e.g. by Giocoli et al. 2013). 

The filtered variance’s strong dependence on non-linear 
structure leads one to think that the skewness of the filtered 
convergence map might be a good discriminator between 
models. The skewness is plotted in Eig. 9 as a function of 
for ah the cosmological models, and once again for three 
source redshifts. As in the case of the variance, the measure¬ 
ments in the EXP003 and EXP008e3 cosmologies are larger 
than those in ACDM, but with the relative residuals that 
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Figure 9. As Fig. 8, but for skewness of the convergence field. 


are a factor of two larger than for the variance. In addition, 
for large angles the relative residuals between ACDM and 
SUGRA003 are very small but increase for small a sig¬ 
nature of the higher small scale clustering present in this 
case, the high concentration of haloes and their dumpiness. 
For a fixed angular scale the difference becomes larger for 
smaller source redshift owing to the build up of non-linear 
structure at late times in the SUGRA003 model. These dif¬ 
ference between models are more evident here than for the 
variance. 

In Fig. 10 we show the variance (left) and the skewness 
(right) of the convergence field as a function of the source 
redshift while fixing the scale of the filter to 6 arcmin which 
is comparable to the typical scale of the central region of 
galaxy clusters at intermediate redshifts. The bottom panels 
again show the residuals of the measurements with respect to 
AGDM. Both the variance and the skewness grow as a func¬ 
tion of the source redshifts because of the increased path 
lengths and the additional structure along the paths. For 
Zs > 0.5 all the models tend to present a constant bias with 
respect to AGDM: positive in EXP003 and EXP008e3 while 
almost vanishing in the bouncing model SUGRA003. Note 
that the relative difference in the skewness is almost double 
than in the variance. While the behaviour of the bouncing 
model SUGRA003 is within the rms of the AGDM measure¬ 
ments, the EXP003 and EXP008e3 are quite distinct with 
a positive bias of about 60 percent and 20 percent for the 
variance and 120 percent and 50 percent for the skewness 
for Zs > 0.5. 

The behaviour of the variance and skewness of the fil¬ 
tered convergence field highlights the possibility of using 
sources at different redshifts and with different smoothing 
scales to investigate the coupling between the dark compo¬ 
nents of the universe. 


4.3 Probability distribution function from source 
redshift distributions 

In previous sections we have discussed the results for the 
variance and the skewness of the convergence field for 
sources located at different fixed redshifts. In this last sec¬ 
tion we explore the one-point distribution function of con¬ 
vergence, shear and magnification given instead a source red¬ 


shift distribution. We do this in order to understand whether 
the simple one-point statistics of the leasing field, and which 
source redshift distribution can help us to understand the 
non-Gaussian nature of the fields in the different cosmolog¬ 
ical models. In our analysis we consider two different source 
redshift distributions that typically correspond to a ground- 
and a space-based weak leasing survey. Eor the ground-based 
case we adopt the source density as a function of redshift as 
computed by Kilbinger et al. (2013) - the data points are 
publicly available on the GEHTLenS webpage^ - while for 
the space-based we consider the parametrization adopted 
by Boldrin et al. (2012) which has been extracted from a 
Euclid-like observation of the Hubble Ultra Deep Eield per¬ 
formed with the code SkyLens (Meneghetti et al. 2008). In 
the top left panel of Eig. 11 we show the two adopted source 
redshift distribution normalizing the GEHTLS one to have 
a total density of sources of about 17 galaxies per arcmin^. 
What is most evident about the two is that the high red¬ 
shift tail in the space-based observation extends to higher 
redshifts while in the ground-based case there are almost no 
galaxies above z = 2. In order to extract the convergence and 
shear catalogues from one light-cone realization we proceed 
in the following way: (i) given a source redshift distribution 
we compute the number of expected sources in a redshift 
interval dz that corresponds to the difference in comoving 
distance between two source planes; (ii) for each dz we ran¬ 
domly assign to each source both a redshift and an angular 
position in the field of view. We do not take into account 
any clustering of the sources which should no be important 
for our purpose. The corresponding lensing quantity (con¬ 
vergence or shear) is then linearly interpolated in redshift 
considering the values computed between two consecutive 
planes at the corresponding angular position of the source. 
The other panels of Eig. 11 show the convergence, shear and 
weak lensing magnification for the different AGDM light- 
cone realizations, each randomly sampled twenty-hve times 
with the according n{zs). Erom the hgure we notice that be¬ 
cause the two adopted source redshift distributions have dif¬ 
ferent shapes, mainly in the intermediate and high redshift 


^ http://www.cfhtlens.org/astronomers/cosmological-data- 
products 
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Figure 10. Variance (left panel) and skewness (right panel) of the convergence field smoothed with a compensated aperture filter with 
6 = Q arcmin as a function of the source redshifts. The bottom panels show the relative residuals of the variance and the skewness with 
respect to the ACDM measurements. The shaded grey region represents the rms of the moments computed on the different realizations 
of the ACDM light-cones. 


tail, the one-point lensing statistics from the same simulated 
light-cones tend to be different. While the shape of the con¬ 
vergence and the magnification are broadened going from 
ground- to space-based observations, the shear distribution 
is shifted toward larger values. 

In Fig. 12 we show the expected convergence, shear and 
magnification probability distributions assuming a space 
based source redshift distribution for the four cosmologies, 
rescaled with respect to the ACDM one. In comparison to 
Pace et al. (2015) we notice that our distributions are realis¬ 
tically more extended in the high-value tails because of the 
different way adopted in doing the ray-tracing, and of the 
larger resolution with which the maps are resolved. While 
in Fig. 6 by Pace et al. (2015) the one-point statistic man¬ 
ifests itself mainly in the different initial power spectrum 
normalisation of the coupled models, in Fig. 12 it is pos- 
sibile to observe also a more pronounced distinction in the 
high value tails. Particularly, in both the PDF of the conver¬ 
gence and the magnification, the small scale clustering and 
the high halo concentration in SUGRA003 model raise the 
high value distribution tails. This is also evident in the shear 
distributions where both SUGRA003 and EXP003 are well 
outside the rms of the ACDM distribution for 7 > 0.075. The 
situation is different for the corresponding PDF extracted 
from the same light-cones considering a ground-based source 
redshift distribution, not shown here. First we notice that 
the source densities in the intermediate and the high red¬ 
shift bins reduce to 11.2 and 0.8 arcmin”^, respectively, 
and that the SUGRA003 model is enclosed within the rms 
of the ACDM for large values. In this case only EXP003 
appears distinguishable from ACDM. The signature of en¬ 
hanced high redshift growth in the considered models is 
much less distinguishable with the ground-based lensing sur¬ 


vey than it is with a space based survey, as also discussed 
in Beynon et al. (2012). 

In Eig. 13 we display the variance of the convergence 
field computed from different light-cone realisations in three 
redshift bins for a space (left) and ground-based (right), 
assuming a CEHTLS source redshift distribution. Black, 
blue, orange and red colors refer to the ACDM, EXP003, 
EXP008e3 and SUGRA003 models, respectively. In each 
panel we also show with a green solid line the noise level 
in the three redshift bins related to the intrinsic elliptic- 
ity distribution Ge = 0.25 and to the corresponding source 
number density Ug contained in the aperture filter (Schnei¬ 
der et al. 1998; van Waerbeke 2000). This is the noise for 
a 25 square degree field. Neglecting systematic errors, the 
noise for a larger survey will go down by roughly a factor of 
one over the square root of the survey area. Erom this plot 
it can be seen that a ground based survey will have diffi¬ 
culty distinguishing between the models, but that a survey 
like Euclid would be expected to clearly distinguish between 
them. 


5 SUMMARY AND CONCLUSIONS 

Several studies have been conducted to understand which 
observables are most suitable to investigate coupled Dark 
Energy cosmologies (Beynon et al. 2012; Giocoli et al. 2013; 
Carbone et al. 2013; Pace et al. 2015). Our aim here was to 
extend the ray-tracing analyses of Carbone et al. (2013) and 
Pace et al. (2015) by investigating whether a tomographic 
slicing of the background sources within high resolution sim¬ 
ulated light-cones might provide additional information to 
distinguish coupled dark energy models from the standard 
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Figure 11. Top left: source redshift distribution expected for space- and ground-based observations. Top right, bottom left and bottom 
right PDF of convergence, shear and magnification extracted from the twenty-five ACDM light-cones randomly sampled eight times with 
the corresponding source redshift distribution. The various color and line style histograms refer to the different considered cosmologies 
(see Fig. 4), the shaded grey region encloses the rms of the measurement performed in the different light-cone realizations in the ACDM 
model. 



Figure 13. Variance of the convergence field computed within the 
different light-cone realisations and cosmologies adopting a space 
and a ground-based (CFHTLS) source redshift distribution. Dif¬ 
ferent data points and colors are the same as in Fig. 10. The solid 
(green) line represents the corresponding noise level associated to 
the intrinsic galaxy ellipticity distribution and the source number 
density. This noise level is for a 25 square degree field. 

ACDM cosmology and possibly the different cDE models 
from each other. Our main results can be summarized as 
following: 

• While the convergence power spectra in the EXP003 
and EXPOOSeS models present only a higher normalisation 


with respect to the ACDM model, the power spectrum in 
the SUGRA003 model exhibits a more complex behaviour. 
Compared to ACDM, in this model structure formation pro¬ 
gresses rapidly at high redshift and slows down at low red¬ 
shift, resulting in a comparable value of erg at z = 0 once 
the normalisation is fixed using the CMB. The large scale 
behaviour of the convergence power spectrum reflects the 
slower growth at low redshifts and it is weakly suppressed. 
This suppression is due to the change of sign of the drag 
term in this class of cDE models, which are characterised 
by a “bounce" of the Dark Energy scalar field whose motion 
changes direction (see Baldi 2012b, for a detailed discussion 
of Bouncing Coupled Dark Energy). In the nonlinear regime 
this effect determines a faster collapse of bound structures 
that results in a higher average concentration of halos and in 
a larger abundance of substructures (see e.g. the discussion 
in Giocoli et al. 2013). This is reflected in an enhancement 
of the convergence power spectrum at small angular scales 
- large 1. 

• The coupling between the dark components can also be 
seen in the lensing signals as a function of source redshifts. 
In particular, if ag is measured from the convergence power 
spectrum, assuming ACDM and fixing Qrn in the fitting, the 
result will change as a function of source redshift in cDE 
cosmologies, i.e. they will be inconsistent with ACDM. 

• The cross-correlation between the convergence for dif- 
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Figure 12. PDFs of convergence (left column), shear (central column) and magnification (right column) extracted from the light-cones 
of the different cosmologies rescaled with respect to the ACDM one, assuming a space based source redshift distribution. Rows refer to 
different intervals for the source redshifts: Zs < 0.5, 0.5 < Zg < 1.4 and 2:5 > 1.4, from top to bottom. Each light-cone has been randomly 
sampled eight time with the according source redshift distribution. 


ferent source redshifts is significantly enhanced in the 
SUGRA003 model while it is slightly suppressed in the 
EXP003 and the EXP008e3 models relative to what is ex¬ 
pected in ACDM, this because of the high concentration and 
the high small scale clustering that manifest as consequence 
of the high and low redshift enhanced and suppressed growth 
rate that characterize the SUGRA003 model. 

• The aperture mass statistic also exhibits signs of the 
coupling between dark matter and dark energy. In particular 
for ^ > 5 arcmin and sources Zg > 1 EXP003 and EXP008e3 
differ by about 50 percent and 20 percent, respectively, as 
compared to ACDM. The skewness of a compensated aper¬ 
ture, measuring the non-Gaussian nature of the convergence 
field, also reflects differences between the various models, 
specifically the higher level of dark matter clustering in the 
SUGRA003 model is evident for small filtering scales. 

• The different models analyzed in the work also result 
in distinct PDE of the lensing signals - convergence, shear 
and magnification - extracted from a space based source 
redshift distribution. The various models not only manifest 
differences in the intermediate values, as discussed by Pace 
et al. (2015), but also in the high value tails. 


Therefore, it emerges from the analyses performed in 
this work that the differences between the coupled dark 
matter-dark energy models and the standard CDM can be 
explained not only in terms of the different normalizations 
of the linear matter power spectrum but also in terms of the 
distinct halo properties and small scale clustering that char¬ 
acterize the dark matter component. These manifest them¬ 
selves in different weak-lensing observables at small angular 
scales and also in the one point statistics taken in different 
redshift bins, given a space based source redshift distribu¬ 
tion that extends to high redshifts. In particular, in this last 
case, both SUGRA003 and EXP003 models appear to be 
easily distinguishable from AGDM by future weak lensing 
surveys. 
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APPENDIX A: LARGE SCALE STRUCTURE 
NOISE IN CLUSTER SHEAR PROFILES 

The results obtained in this paper can also be used to evalu¬ 
ate how the level of noise produced by the large scale struc¬ 
ture on the estimates of galaxy cluster masses from weak 
lensing data (Bahe et al. 2012; Hoekstra et al. 2012; Merten 
et al. 2014; Giocoli et al. 2014; Sereno et al. 2014; von der 
Linden et al. 2014,b) are dependent on the cosmological 
background. 

Several algorithms have been developed to invert the 
lens equation and derive the cluster’s projected mass (Julio 
et al. 2007; Merten et al. 2009; Zitrin et al. 2011). Since 
galaxy clusters have a matter distribution that on average 
can be described by a well defined density profile (Navarro 
et al. 1996, 2004), a direct way to weight them is to fit the 
(spherically binned) measured shear profile adopting the cor- 
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responding theoretical prediction (Bartelmann 1996; Giocoli 
et ah 2012a). 

However, the light traveling from sources located be¬ 
hind a galaxy cluster is deflected not only by cluster matter 
distribution, but also by all the matter it encounters along 
its trajectory. This uncorrelated matter density distribution 
contributes to the lensing signal, and affects the shear mea¬ 
surements (Schneider et al. 1998; Hoekstra 2003; Hoekstra 
et al. 2011). 

This lensing noise, which depends on the source red- 
shift and on the angular size 0 of the annulus we are look¬ 
ing at, can be computed analytically from the convergence 
power spectrum adopting the formalism developed by Hoek¬ 
stra (2003): 

olssio) = ^ IdiiPGiuiiW), (Al) 

where J 2 represents the second-order Bessel function and is 
a particular choice of the aperture mass statistic averaged 
over an annulus ranging from {0 — 50/2) to {0 -\- 50/2). 

In Fig. Al we show, for the four considered cosmolo¬ 
gies, the noise produced by the LSS on spherically averaged 
shear profiles. The results are presented for three typical 
scales (^=1, 10 and 20 arcmin) and considering sources up 
to redshift Zs — A. As already done in the case of the aper¬ 
ture mass statistics, the results are directly extracted from 
the convergence maps that have been smoothed using the 
appropriate J 2 kernel. 

From the figure we notice that the noise produced by 
large scale structures is higher in the cDE models with re¬ 
spect to the ACDM cosmology. This is particularly true 
for the EXP003 and EXP008e3 models: for ^ > 10 arcmin 
and for sources with > 0.5 they present values that are 
about 20 and 10 percent larger than the one measured in 
the ACDM model, respectively. Eor the SUGRA003 model, 
a smaller increase of noise is evident only for 0—1 arcmin 
and/or for sources at low redshift. We underline also that 
the rms of the measurement ACDM decreases as a function 
of redshifts because of a combination of (z) the choice filter, 
{ii) the field of view and (iii) the clustering of the haloes 
as a function of redshifts - at high redshift the universe is 
more homogenius. 
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Figure Al. Large scale structure noise computed adopting a particular choice of the aperture mass statistic, averaged over an annulus 
ranging from (6 — 66/2) to (6-\-S6/2)^ as a function of the source redshift measured from the convergence maps extracted of the light-cones 
of the four considered CoDECS simulations. Left, central and right panel refer to different angular scales 6] colors and line styles are as 
in Fig. 8. 
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